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ABSTRACT 

Context. Results from the theory of radiatively driven winds are nowadays incorporated in stellar evolutionary and population syn- 
thesis models, and are used in our interpretation of the observations of the deep Universe. Yet, the theory has been confirmed only 
until Small Magellanic Cloud (SMC) metallicities. Observations and analyses of O-stars at lower metallicities are difficult, but much 
needed to prove the theory. 

Aims. We have observed GHV-62024, an 06.5 Illf star in the low-metallicity galaxy IC 1613 (Z ~ 0.15Z o ) to study its evolution and 
wind. According to a previous preliminary analysis that was subject to significant restrictions this star could challenge the radiatively 
driven wind theory at low metallicities. Here we present a complete analysis of this star. 

Methods. Our observations were obtained with VIMOS at VLT, at R » 2000 and covered approximately between 4000 and 7000 A. 
The observations were analysed using the latest version of the model atmosphere code FASTWIND, which includes the possibility of 
calculating the N in spectrum 

Results. We obtain the stellar parameters and conclude that the star follows the average wind momentum-luminosity relationship 
(WLR) expected for its metallicity, but with a high value for the exponent of the wind velocity law, /?. Comparing this with values of 
other stars in the literature, we suggest that this high value may be reached because GHV-62024 could be a fast rotator seen at a low 
inclination angle. We also suggest that this could favour the appearance of the spectral "f '-characterictics. While the derived [3 value 
does not change by adopting a lower wind terminal velocity, we show that a wrong V„ has a clear impact on the position of the star 
in the WLR diagram. The N and He abundances are very high, consistent with strong CNO mixing that could have been caused by 
the fast rotation, although we cannot discard a different origin with present data. Stellar evolutionary model predictions are consistent 
with the star being still a fast rotator. We find again the well-known mass-discrepancy for this star. 

Conclusions. We conclude that the star follows the WLR expected for its metallicity. The results are consistent with GHV-62024 
being a fast rotator seen close to pole-on, strongly contaminated at the surface with CNO products and with a wind structure altered 
by the fast rotation but without modifying the global WLR. We suggest that this could be a general property of fast rotators. 

Key words. Galaxies: individual: IC1613 - Stars: early-type - Stars: fundamental parameters - Stars: mass-loss - Stars: rotation - 
Stars: evolution 



q ! 1. Introduction iKudritzki & Puis! d2000t) and lPuls. Vink & Naiarro I d2008l) . One 

(~0 of the most important predictions of the theory is that the modi- 

^ . Three main parameters determine the structure and evolution of fied wind momentum rate of the stellar wind (MWM? the prod . 

L" ! a star: lts mass ' an 8 ular momentum and metallicity. For massive uct of mass . loss rate> wind termina l velocity and square root of 

m ;> , stars the first two parameters change rapidly during their lives the stellar radius> which we denote as = M Vca R^ 2 follow- 

k> ■ because of stellar winds, whose strength is determined by the ing the Ufje in the literature) should depend almost solely on the 

; third parameter. A reliable theory of stellar winds for massive stellar i uminosity and metallicity. The VLT-FLAMES survey of 

Ctf , stars is accordingly crucial for our predictions about how these masfjive starfj (E vansetalj|2005, P.I. S. Smartt) was a project 

stars evolve and therefore for our predictions of their properties aimed prec iselyat testing this prediction, which was nicely con- 

(lomzing flux, mechanical energy delivered, chemical composi- firmed by Mokiem et al. (2007b) by analysing O-type stars in 

tion, etc.) at a given moment in their lives and final fate (type of the MOky Way> the LMC and me S MC. However, the second 

supernova and remnant left behind). part of the pro j ect; the comparison of the predictions of stellar 

The theory of radiatively driven winds, used to explain the evolutionary models including mass-loss and rotation with the 

observed winds in massive stars and whose predictions are in- ste n ar properties revealed problems when comparing the nitro- 

corporated in stellar evolutionary and population synthesis mod- gen abundances (iHunter et al. II2008I) . indicating that other phys- 

els, is based on the work by jCastor, Abbott & Klein | (|1975|), ical proce sses may be at work. This started the VLT Tarantula 

with improvements by iPauldrach, Puis & Kudritzki| dl986|) and FLAMES Survey (VFTS), a detailed investigation of the mas- 

|Friend& Abbott | dl986|). Recent reviews can be found in s i ve s tellar population of 30 Doradus in the LMC currently under 

way dEvans et al. II201 ll) . There a deeper study of the behaviour 

Send offprint requests to: A. Herrero (e-mail: ahd@iac . es) of the LMC O-stars winds will be possible, but without extend- 

* Based on observations obtained at the ESO VLT for Programmes m „ me aV ailable metallicity baseline 
078.D-0767 and 080.D-0423. 
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Very many O-stars have been analysed in the Magellanic 
Clouds (see, among others, pouret et al. 1 12.003b iHillier et al. I 
20031: IMokiem et al. Il2006l IMokiem et al. Il2007t iMassev et al. I 
2004 or IMassev et al. 1120 05?). However, going to lower metal- 



licities implies a large distance jump. Yet, many star-forming 
galaxies in the Local Universe and presumably most stars in the 
early epochs of the Universe have metallicities below that of the 
Small Magellanic Cloud. It is therefore natural to try to go one 
step further and analyse stars at even lower metallicities, study- 
ing their wind behaviour and stellar evolution. Although going to 
galaxies farther away and with lower metallicity than the SMC 
will result in an increasing difficulty to obtain good quality spec- 
tra that can be analysed, present day telescopes and instruments 
clearly allow for the next step. At metallicities of about half that 
of the SMC and distances of one to a few megaparsecs it is 
still possible to perform reliable analyses. This has been done 
using the visua lly brightest objects, the A- and B- supergiants, 
bvlUrbaneia et al. I d2008l) in the WLM (0.96 Mpc). [Evans et al. I 
( 2007b in NG C3109 (1.3 Mpc ) ICastro et all (l2012h in NGC55 
(2.0 Mpc) or [Urbaneia et al. I (|2005|) in NGC 300 (2.2 Mpc). 
iBresolin et al. " 3200 ll) have even analysed individual A super- 
giants in NGC 3621 at 6.7 Mpc. The hotter O-stars, however, are 
visually fainter and much more difficult to observe at far dis- 
tances. For this reason, there have been no spectroscopic analy- 
ses of individual O-stars located in low-Z galaxies using modern 
standard techniques (i.e., analysing optical spectra by means of 
NLTE, spherical models wi th mass-loss), not even in the Local 
Group , until very recently dHerrero et al. II20TTI iTramper et al. I 
I20T1 . 

IC 1613 is a dwarf irregular galaxy in the Lo cal Group with 
a dist ance modulus of (m - M)o = 24.31 + 0.06 ( Dol phin et al. I 
120011) . It has a metallicity of log( Q/i/) + 12 = 7.80 ±0 .10 as de- 
termined from its B-supergiants (Bres olin et al. 20071). wherea s 
nebular studies vary between 7.60 and 7.90 dLee et al. 2003). 
The metallicity of the B-supergiants indicates Zic 1613 — .13 ± 
0.07 Z (using solar abundances from Aspl und et al. 12 009). This 
value is below the value deri ved by most autho rs for the SMC 
O- and B-stars. For example. [Heap et ai~l (|2006[) obtain Zsm c = 
0.2 Z from an analysis of O-stars and lDufton et al. I (120051) ob- 
tain Zsmc = 0.25 Z from an analysis of B-supergiants. Th e 
situation is similar for other authors, e.g. iBouret et al. I (120031) . 
who obtained Zsmc = 0.2 Z . Therefore the stellar metallicity in 
IC 1613 is below the SMC value and the study of its stars repre- 
sents a step forward in our analyses at low metallicities. 

IC 1613 shows a recent and intense burst of massive star for- 
mation, particularly in its NE part. We have recen tly published 
a new catalogue of OB associations in IC 1613 dGarcia et al. I 
2009, identified as [GHV2009]) and their physical properties 
dGarcia et al. 2010) as part of our effort to carry out an in-depth 
study of the young population of IC 1613. 

We also obtained spectra of some stars in this galaxy us- 
ing VIMOS at VLT. In the field of IC 1613 we found the star 
62024 in the GHV2009 catalogue (hereafter we identify the star 
as GHV-62024), an Of star. This star is ideally suited to study 
the stellar evolution and winds in IC 1613, because it displays 
strong Nm lines and a P-Cyg ni profile in Hen 4686 . A prelim- 
inary analysis of this star by iHerrero et al. I (1201 lb resulted in 
a too high logD mom as expected from its metallicity and lumi- 
nosity, which would challenge the theory of radiatively driven 
winds at low metallicitie s. An even mo r e serio us challenge is 
the very recent work by ITramper et al. I (1201 lb . These authors 
analysed six O-stars in low-Z galaxies (four in IC 1613 and one 
in WLM and NGC 3109). The four stars in IC 1613 tend to show 
too high values of D mom , although they are not inconsistent with 



the theory because of the large error bars. But the two other 
stars have too high D mom values, not consistent with theoreti- 
cal predictions. However, although these two recent works seem 
to point towards a problem between observations and theory, we 
note that the analyses are subject to strong assumptions. 

In the present work, we carry out a detailed analysis of GHV- 
62024, showing that when we drop the assumptions made by 
IHerrero et al. I ( 20TT1) . we obtain a MWM value consistent with 
the theoretical expectations, although the analysis opens alterna- 
tive questions. The observations and data reduction are presented 
in Sect. [2] The spectrum is described in Sect. [3] The spectral 
analysis follows in Sect. [4] and a discussion of the stellar wind 
momentum rate is given in Sect. [5] The evolutionary status and 
other properties are discussed in Sect. |6] Finally, the conclusions 
are presented in Sect. [7] 



2. Observations and data reduction 



Details of the observations are given in IHerrero et al. I d20l0h . 
We repeat here the main points. Observations were performed 
with VIMOS at VLT in MOS (multi-object spectroscopy) mode. 
The HR-blue and HR-orange gratings were used, resulting in a 
resolution of R ~ 2050 and R ~ 2150, respectively, and a spec- 
tral wavelength coverage roughly from 3870 to 7240 A (de- 
pending on the star position). We had 19 observing blocks for 
the blue part of the spectrum, each with three individual expo- 
sures, and ten observing blocks for the red part, each with two 
individual exposures. Thus, in total 3x19 blue spectra were ob- 
served from October 5 to November 5, 2007, and 2 x 10 red spec- 
tra from November 5 to November 13, 2007. Integration times 
for individual exposures were 650 sec for the blue spectra and 
1005 to 1110 sec for the red spectra. The three (two) consecu- 
tive exposures in the blue (red) were coadded. The coadded spec- 
tra were extracted with standard IRAFQ procedures after wave- 
length calibration and cosmic ray removal. Next, a barycentric 
correction and a co rrection for the systemic velocity of IC 1613 
(-234 + 1 km s -1 , iLu et al. Ifl993h were applied. Finally, after 
checking that no significant radial velocity variations are present, 
all 19 (10) resulting spectra were coadded to a single blue (red) 
spectrum. The resulting spectra, with an S/N of 85 both in the 
blue and red wavelength regions, were then rectified. 

Figure Q] shows the finding chart of GHV-62024. The star 
is located at a(2000)= Olh 05m 00.72s, £( 2000)= 2° 8' 49". 1 
(from the astrometry by Garcia et al. 1 120091) . It clearly lies in a 
zone of intense recent star formation, with strong ionized bub- 
bles visible in Fig.Q] Table Q] gi ves the optica l photo metric data 
of the star from the catalogue of Garci a et al. 1 (120091) . 



Table 1. Photometric data for GHV-62024 from iGarcia et al.l 
(2009). 
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R I 


18.69 


19.44 


19.60 


19.76 19.68 



1 IRAF is distributed by the National Optical Astronomy 
Observatory, which is operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with the 
National Science Foundation. 
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INT-WFC, R-band-; 



^LT-VIMOS, R-band 



62024 N_ 



Fig. 1. Finding chart for GHV-62024. The left part is an image 
obtained with the WFC at the INT in the R band. The region 
imaged (5x5 arcmin) contains the IC 1613 intense star-forming 
region, with many ionized bubbles. The right-hand side is a VLT- 
VIMOS image in the same band, at a larger scale (indicated by 
the scale bar at the right edge). The orientation is the same in 
both images. 



3. Description of the spectrum 

Figures [2] and [3] show the observed coadded blue and red spec- 
tra of GHV-62024, respectively (together with the best-fit model 
that we describe below). The stellar spectrum is dominated by 
the H and He lines, plus a conspicuous N in emission at 4630-40 
(but without N iv) and a weak Si iv 4089 (but Si iv 41 16 could not 
be detected). The core of the Balmer lines show an intense neb- 
ular contamination, but fortunately the wings are well defined 
(except for H a ). The red spectrum has a poorer quality, and no 
individual features could be unambigously identified, apart from 
H and He lines. 

The star has a Hen 4541 line stronger than Hei 4471. 
By measuring the equivalent widths of the se lines (log(W) = 
log(£W447i) - log(£W4542) = -0.18, see IConti& Alschulerl 
ll97lHMarkovaet al.l l20in. and comparing with standards we 
classify the star as 06.5. Because of its Nm 4634-40 and Hen 
4686 emissions we assign it the f-class. The luminosity criteria 
of O stars are based on metal lines that are very faint in our spec- 
trum and are usually defined at other metallicities. We assign 
the star a luminosity class based on the strength of the Balmer 
lines. Thus the star is classified as 06.5 Illf. The red spectrum 
is heavily contaminated by nebular and sky emission and conse- 
quently H ff yields no information on mass-loss rate. Therefore, 
the P-Cygni profile at Hen 4686 is a very important signature. 
The 19 individual spectra have an average SNR of 20 (with a 
large scatter). We then inspected the spectra of individual ob- 
serving blocks to ensure that the P-Cygni profile is not caused 
by unidentified cosmic rays, spikes, or just random noise vari- 
ations. We also inspected the nebular emission around the star, 
and confirmed that there is no emission in He n 4686. We con- 
clude that the P-Cygni profile is a stellar feature. 



4. Spectral analysis 

We analysed GHV-62024 with the newest FASTWIND 
version, which inc l udes dielectronic recombination 
dRivero-Gonzalez et al. I 1201 lal) . The analysis was based 
mainly on U 6 , H y and H^, Hei 4471, Hen 4541, 4686 and the 
N in emission lines at 4634, 4640 A. Although all line profiles 
are sensitive to all parameters, the main diagnostics for gravity 
was the wings of the low Balmer lines (except H„ which was 
not used). For the temperature we mainly used the ratio He i 





t 







Fig. 3. Same as Fig. [2] red spectrum. 



4471 to Hen 4541, giving some weight also to the Hei 4387, 
4922 lines (which point to a temperature slightly cooler than 
Hei 4471). The ratio of Nm 4634-40 to Niv 4058 also was a 
secondary indicator (the N iv line is not present, thus giving an 
upper limit for the temperature). Without a useful H a profile, 
the mass-loss rate and the exponent of the velocity law, /3, were 
determined using the Hen 4686 line. This is crucial, because 
there is a degeneracy between the B exponent and the derived 
mass-loss rate. In Fig. [4] we see the fit to the Hen 4686 profile 
using the value of B= 2.0 and B- 1.0. From the figure it is 
clear that the higher value fits the observed profile better. It 
turned out to be impossible to find a good fit to the red wing of 
He ii 4686 with a low B value. The abundances of N and He, 
given in Table |2l were fixed by fitting the profiles of all lines. 
The first one is given in the usual logarithmic scale in which 
the H abundance is 12, while the He abundance is given by 
number, by Fee = (Nhs/Nh). The red spectrum was not used for 
the analysis, but we note that it is consistent with the adopted 
parameters (see Fig. [3). 

The final stellar parameters are given in Table |2]and the best 
fits are shown in Figs. [2] and [3] For the analysis we adopted a 
projected rotational velocity of 80 km s and a radial-tangential 
macroturbulent velocity of 50 km s -1 . Although we were un- 
able to determine these values from the line profiles, at a spec- 
tral resolution of R « 2000 they have little impact on the pa- 
rameter determination, because the line broadening is limited 
by the spectral resolution or the Stark broadening. The adopted 
metallicity for the model calculations was Z = 0.15 Z Q at the 
beginning and we changed to Z = 0.13 Z for the models 
close to t he final stellar param eters (adopting the metallicity ob- 
tained by Bres olin et al. Il2007l for B-supergiants, although this 
difference in Z produced no appreciable change in the spec- 
trum). Unfortunately, there are no features in our spectrum that 
allow us an independent metallicity determination. The radius 
was determined by fitting the visual magnitude by integrating 
the flux of our final model in the Johnson V-filter and after 
correcting for the d ist ance to IC1613 and the reddening (see 
iHerrero et al. 1119921 orlRepolust et al. 112004b . We used the pho- 
tometry from iGarcia et al. I (|2009), quoted in Tab. Q] The final 
model has an intrinsic (B - V)o = -0.26, while the observed 
one is (B — V) = -0.16. Assuming a standard extinction law 
with R v = 3.1 and with the adopted distance modulus of 24.31 
we finally obtained an absolute magnitude M v - -5.03, which 
results in a radius R — 1 1 . IR & and using the derived effective 
temperature, in a luminosity log(L/L Q ) = 5.29. These values 
are consistent with a slightly underluminous 06.5 giant, which 
reduces the possibility that the star is actually a binary system 
with two comparable components that could affect the spectral 
analysis. 
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Fig. 2. Observed blue spectrum of GHV-62024 (red) and the adopted best fit (black). The main He and N lines used in the analysis 
have been marked. Stellar parameters are given in Table |2] The emission in the red wing of H r is a sky line. 



We have little information on the terminal velocity. We ini- 
tially adopted Voo/ Vesc = 3. 4 and V^iZ) oc Z° 12 as given by 
ICastro et al.l d2012l) (see also ICastro 11201 Ol) , who used data from 
the literature (the m etallicity dependence has been taken from 
iLeitherer et al. Ifl992l) . We note that this ratio is higher than the 
more commonly used one of 2.6, which would result in a termi- 
nal velocity 30% lower. This has little effect on the analysis, but 
may be important for the MWM (see next section). We adopted 
a wind terminal velocity of 1800 km s _1 . This is consistent with 
the absorption part of the He n 4686 P-Cygni profile, which in- 
dicates a lower limit of »T200 km s _1 , although a more precise 
determination of Voa from this line is not possible. The mass-loss 
rate is derived mainly from the emission profile of He n 4686. 
We obtain M = 4.3 x 10~ 7 M Q yr" 1 . 

We estimated the error in r e ff to be ±1500 K and in logg ± 
0.15 dex, which translates into an uncertainty of +0.6 in R/R Q 
and +0.09 in log(L/L ). The spectroscopic stellar mass (the 
mass obtained using the values of log g and radius obtained from 
the analysis of the spectrum, and corrected from the centrifugal 
force) results in M sp = 18. 3*ZJ M . The error in the helium 
abundance is estimated to be ±0.04 and in log(N/H) ± 0.30 dex. 
The uncertainty in V x is large, ±600 kms -1 . The uncertainties in 
P and M are more difficult to estimate because the upper limit 
of jS (and thus the lower limit of M) is not well constrained. 
However, if we consider p fixed, the uncertainty in M (which 
is then a lower limit) is ±0.30 dex. 

Table 3. Synthetic photometry for the model given in Table [2] 
Values are absolute magnitudes. 



u 


B 


V 


R I 


-6.37 


-5.29 


-5.03 


-4.84 -4.61 



5. The wind momentum of GHV-62024 

The first remarkable result of the previous section is the high 
value we obtain for /3, the exponent of the wind velocity law. 




WAVELENGTH (A) 



Fig. 4. Fit to the Hen 4686 profile using different values for 
the terminal velocity and the p exponent of the wind velocity 
law. The mass-loss rate was also varied in each case to ob- 
tain a reasonable fit while keeping all other parameters fixed. 
Red, dotted line: = 2.0, Voo = 1800 kms 4 ,M = 4.3 x 
10~ 7 M yr~' (adopted model); blue, dashed line: p = 1.0, - 
1800 km s _1 ,M = 1.1 x 10~ 6 M yr _1 ; orange, dot-dashed line: 
= 1.0, Voo = 900 kms"',M = 5.5 x 10~ 7 M yr 1 . Black: 
Observed spectrum. 



The value obtained (J3 = 2.0) is much high er than usual values 
for early and m id Of stars (0-7-1 .0, se e e.g. lRepolust et al.ll2004l 
for the MW or lMokiem etaTll2006l for the SMC). However, p 
and M have different effects on the line profiles. A higher p 
means a slower acceleration in the inner wind (where the He n 
profile is formed) and therefore, because of the continuity equa- 
tion (M = Ajir 2 p{r)v(r)) for a fixed mass-loss rate and terminal 
wind velocity, at a given point r the velocity is lower and the 
density higher than for a smaller p (see Fig. |5j. This implies 
more emission at low velocities, but also a lower velocity at a 
given density. The minimum density required to contribute to 
the emission profile is reached at lower velocities in the high p 
model. As a consequence, the line profile in this model has more 
emission at low velocities compared to the high velocities. The 
result is a steeper red wing, fitted with a lower mass-loss rate in 
the high p model. 
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Table 2. Stellar parameters obtained here for GHV-62024 with FASTWIND. T e ff is given in K, g in cm s~ 2 , M in solar masses per 
year and V M in km s 



Star ID SpType 7~ eff logg R/R Q log(L/L ) M sp /M M f3 Y Hc log(N/H) + 12 

62024 06.5 Illf 36500 160 TTl 5T29 183 4.30x10-' 2l) 1800 018 830 



E - 



z 

LU 
Q 



-0.2 



0.0 



0.2 



0.4 

V/Vinf 



0.6 



0.8 



Fig. 5. Run of the density versus wind velocity for three models. 
The first one (solid line) has M = 1.0 x 1(T 6 M Q yr-\/3 = 1.0 
and the second one (dashed line) has M = 4.3 x 10~ 7 M Q yr~ l ,l3= 
2.0. In both cases, the mass-loss rate is the one need to fit the 
emission peak of Hen. The third model (dot-dashed line) has 
M = 4.3 x 10" 7 M e yv~ l ,f3 = 1.0. All other parameters are the 
same as in our final model. 



In a preliminary analysis of GHV-62024, iHerrero et al. I 
(1201 11) assumed a standard value of f3- 0.9 for its wind. As a 
result, the M WM rate of this s tar was much higher than pre- 
dicted by the IVink et al~l (12001 1) relationship for its luminosity 
and metallicity (note that if the star is contaminated by a com- 
panion, which is always a worry at far distances, the actual lumi- 
nosity would be lower and thus the discrepancy with the theory 
larger). As in dicated in the introdu ction, this result is similar to 
that found bv lTramperet al. I (1201 lb for stars in IC 1613, WLM 
and NGC 3109, who also adopted ft (with values from 0.80 to 
0.95, depending on luminosity class). As surprising as the /3 
value obtained in this work may look, it solves the problem of 
the high MWM of GHV-62024, because it reduces the mass-loss 
rate while keeping all other parameters (nearly) fixed. Figure [6] 
illustrates this point. In the figure we plot the present lo g(D morn ) 
value (28.22) and that obtained bv lHerrero et al. 1(1201 lb (28 .69), 



together with the values given by Mokie m et al. ( 2007bl) for 



Galactic and SMC stars. The stellar luminosity has also changed 
because we derive slightly different stellar parameters as a con- 
sequence of the new The new value for the MWM rate fits 
the cannonical relationship defined by the SMC stars for un- 
clumped winds. We note, however, that because the metallicity 
of IC 1613 is slightly lower than that of the SMC, we could ex- 
pect the log(D mom ) of GHV-62024 to be slightly lower as well. 

The high (3 value obtained here is unusual, but there are 
precedents in the literature. The most similar ca se is probably 
that of AzV 83, analysed bv lHillier etaTI (120031) . This is an 07 
If star (as compared to 06.5 Illf for GHV-62024) for which the 
authors also found a velocity law exponent f3= 2.0. We do not 
include clumping in our analysis, but Hillier et al. did, therefore 
this is not the cause of the large This high value of (3 for such 



an early star led the authors to suggest that AzV 83 is a fast rota- 
tor seen pole-on. The peculiar value for /3 would be the result of 
a mixture of photons from the asymmetric wincfl Changing ac- 
celeration conditions along the wind may mimic a high /? value 
(Owocki, private communication). Extending this argument, we 
may expect that the equatorial and polar winds are subject to 
different acceleration conditions, and accordingly that the global 
profile mimics a velocity law with a high p. 

There a r e mo re similarities between the two stars. 
iHillieret all (120031) obtained F H e = 0.20 and 12 + log(N/H) = 
8.41, compared to our F H e = 0.18 and 12 + log(N/H) = 8.30. 
In both cases, the high abundances could be attributed to the 
proposed fast rotation. The main difference between both ob- 
jects lies in the wind terminal velocity. AzV 83 has a modest 
Vco = 940 km s -1 as derived from its UV spectrum. This value is 
much lower than the 1800 km s _I we estimated for GHV-62024 
by scaling its escape velocity. The difference would have been 
smaller had we also used the escape velocity to determine the ter- 
minal wind velocity of AzV83. The escape velocities of the two 
stars are not so different: 436 km s _1 for AzV83 and 666 km s _1 
for GHV-62024 (once more, it is true that the devil is in the de- 
tails). Thus, when estimated from the respective escape veloc- 
ities, the terminal velocities of both stars would differ by only 
-30%. 

Could it be that the high f3 value derived here for GHV-62024 
is the result of a too high estimation for V^? As was shown by 
iPuls et al.l (Il996l) . it is primarily the Q-parameter ( Q = (y ^ )3 / 2 ) 
that determines the net wind emission. Therefore, for a given 
wind emission, if Voo increases M has to increase as well. And 
given the behaviour of density and velocity explained above and 
illustrated in Fig. [5] this could result in too much emission at 
high velocities. Could a lower V M result in a lower /3 value? The 
answer is no (see Fig. |4). Although a given stellar parameter has 
an effect on the whole profile of a given spectral line, its impact 
on different parts of the line is also different. Consequently f3 
primarily determines the steepness of the red wing emission (be- 
cause a higher /3 means a lower wind acceleration, and therefore 
the inner layers are denser and slower, resulting in more emis- 
sion at a given wavelength but extending to lower velocities). 
On the other hand, the main impact of Voo is in the extension of 
the lines profiles, primarily in the absorption part of the P-Cygni 
profile. In Fig.|4]we show an example with different values of Voo 
and j6. We see that we can reduce the red emission by decreas- 
ing the terminal velocity at f3- 1 .0, but then neither the slope of 
the red wing nor the extension of the blue absorption would fit. 
We conclude that our estimation of the wind terminal velocity 
cannot be wrong by more than 30%, which is consistent with the 
lower limit indicated in Sect. [3] and that this would not affect the 
derived /3 significantly. 



2 Hillie r et al. I d2003 | ) adopted the equ atorially compressed wind 
model from Biorkman & Cassinelli ( 1993), i.e., the equatorial wind is 
denser than the polar one, in contrast to current views, which favour a 
a faster and denser wind from the pole and a slower, thinner wind from 
regions closer to the equator. However, this difference is not important 
for the argument presented here. 
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AzV 83 is one of the stars conside red by Mokiem et aTl 
(l2007bl) (who adopted the values given bv lHillier et al. Il2003l af- 
ter correcting for the clumping adopted by these authors) and 
appears in Fig. [6] in a position fully compatible with the rest of 
the stars. Again, a more standard ft value would result in a MWM 
value higher than the average relation. 

However, a too high wind terminal velocity would have a 
clear impact on log(D mom ). As log(D mom ) = log(MVcoR 0S ), 
using the above definition of Q we obtain log(D mom ) = 
2.5^^00 + 2 log/? + log Q. For a fixed pair (Q,R) we obtain 
that a change of 30% in V M results in a change of 0.3 dex for 
log(D mom ). This would bring the point of GHV-62024 in Fig. [6] 
close to the value expected for the IC 1613 relationshi p. This is 
what o ne would expect even for a fast rotator, because M aeder I 
(2001) has shown that the WLR (and in turn the position of the 
star in Fig. [6j is not altered by fast rotation. 

Therefore, we conclude that these two stars are consistent 
with a scenario in which they are fast rotators seen pole-on, their 
wind momentum is what is expected from their luminosities and 
metallicities. We suggest as a working hypothesis for future in- 
vestigations that the fact that their ft values are higher than usual 
may be a consequence of fast rotation. 

Can we also observe higher ft values for known fast rot ators? 
Two fast rotators in the sample of Mo kiem et al. I (1200 7b) lie at 
positions close to GHV-62024 in the W LR diagram. NGC346 - 
010 is an 07 IIIn((f)) star analysed bv IMokiem et al. I J2006). 
NGC 346-010 is in fact quite similar to GHV-62024, except for 
the higher mass-loss rate and He abundance of the latter. Fig. [7] 
offers a comparison of the spectra of NGC346-010 and GHV- 
62024, with the last one convolved with an additional rotational 
profile with a velocity of 200 km s (to simulate the high rota- 
tional velocity of NGC346-010, but taking the much lower spec- 
tral resolution for GHV-62024 into account). In the figure we 
appreciate that both spectra are very similar, except for the emis- 
sion at He n 4686 and N m 4630-40, a higher ratio He n 4542 / 
He i 4471 and the nebular emission in H r . The similarity in the 
rest of the spectral features indicates that the GHV-62024 emis- 
sion in He n 4686 is a real feature. This figure can be c ompa red 
with the spectrum of AzV 83, shown bv lHillier et al. I (12003b in 
their Fig. 2. There we see the strong emission in N 4634^-0 and 
He 4686 (this time the line is fully in emission). Thus, the three 
stars look quite similar, except for these emission features (and 
the small differences in the other He profiles related to effective 
temperature, which could also be affected by differences in at- 
mospheric extension). Because of the lack of UV data, the wind 
terminal velocity of NGC346-010 has been again estimated from 
the es cape velocity (with V M — 1832 km s" 1 in Mokiem et al. I 
2006, later revised to 1486 km s" 1 in lMokiem et al Il2007bl) . A 
standard value, ft = 0.8, was adopted for the exponent of the ve- 
locity law. For this fast rotator, the authors suggest the presence 
of an equatorial disk (because of some difficulties in the line pro- 
file fitting, but without other evidence) although on the whole the 
fit to the H/He optical s pectrum is satisfactory. The second fast 



rotator in the sample of Mo kiem et al. | (|2007bl) is AzV 2 96, a 



fast rotating 07.5 V((f)) analysed by |Massey etaT] (l2004l) . In 
this case, the terminal velocity has been obtained from UV spec- 
tra ( Voo = 2000 km s _1 ). Aft = 0.8 has also been adopted for the 
exponent of the velocity law. In this case the optical fit is poor 
and the authors suggest that the star might be a binary, which 
would also explain the relatively large radius needed to fit the 
observed magnitude and the high radial velocity observed. 

In both cases, a log(D mom ) slightly larger than the cannoni- 
cal relationship is derived. Unfortunately, all lines are strong in 
absorption including H„, and therefore the ft exponent could not 



be accurately constrained. Consequently, the mass loss rate in 
this case may actually be an upper limit. A higher/? value would 
therefore bring these two fast rotators into better agreement with 
the average WLR for their metallicity. 

Therefore, we conclude that the analysis of these stars in- 
dicates that the velocity law may be altered in some stars. 
Anticipating data from the next section (where we discuss the 
N and He abundances), we suggest that fast rotation is a prob- 
able cause for this alteration. Note, however, that the evidence 
is indirect and based on very few objects. Moreover, of three 
stars with known "f" signature at low Z (AzV 83 and NGC 346- 
001 in the SMC and GHV-62024 in IC 1613), two have been 
suggested to be seen pole-on. This is quite improbable, unless 
the f-phenomenon at low Z is favoured by this pole-on view. If 
that were the case, the enhanced polar wind would be respon- 
sible for the presence of the f characteristics. This possibilit y 
agrees well with the findings bv lRivero-Gonzalez et al. I (1201 lal) . 
who found that the f-phenomenon is related to wind strength and 
mi ght therefore de c rease at low metallicity. 

Mok iem et al. I d2007bl) have analysed other fast rotators in 
the SMC, but they give no further clues because the authors 
could only derive upper limits that are above the mean WLR 
for the SMC. However, it is interesting to note that we do 
not find this si t uation in th e Milky Way: the valu es derived by 
IMokiem et aTl (l2007bl) and lRepolust et al. I (120041) . whether ac- 
tual values or upper limits, are close to or below the average 
WLR for the Milky Way and we are not aware of significantly 
high ft values needed to fit the spectra of Galactic O stars, al- 
though they are found for the slower winds of B-supergiants 
(Kudritzk ietal.il 19991; ICrowtheretaT~ll2006l: iMarkova & Puis] 
2008). We note that determining ft from optical spectra of rapidly 
rotating stars seen edge-on is in most cases very difficult because 
of the lack of spectral features sensitive to this parameter. 

We have not taken clumping into account in our analysis 
and therefore t he derived mass-loss rate has to be considered as 
an upper limit (lPuls. Vink & Naiarro 112008b . However, as far as 
we compared ours with other results and theoretical predictions 
that also use unclumped winds, our conclusions remain valid. 
Of course, there might be relative differences in the clumping 
among stars with different properties, in particular the winds of 
fast rotators might be clumped in a different way than those of 
slowly rotating stars. However, studies about the clumping dis- 
tribution are at their very beginning, and therefore we have no 
data yet to discuss this point. 

6. The evolutionary status of GHV-62024 

In Fig. [8] we place GHV-62024 in the Hertzspmng- Russell dia- 
gram and compare it with evolutionary tracks from iBrott et aTTI 
(1201 lh for the SMC. We chose the tracks with an initial rota- 
tional velocity V rot = 410 km s _1 because they have the high- 
est rotational velocity consistent with the star locus. For higher 
rotational velocities, the corresponding tracks evolve bluewards 
already from the ZAMS. 

We see that according to these tracks the star was bom with 
an initial mass of ~30 M Q . This result does not change signifi- 
cantly when using other tracks, even those without initial rota- 
tional velocity. The present evolutionary mass derived from the 
adopted tracks (at V rot; ini = 410 km s~') is 28.6M . A compar- 
ison with the spectroscopic mass reveals then the well-known 
mass-discrepancy. 

The most interesting parameter for the evolutionary status of 
GHV-62024 is its nitrogen abundance. We obtained a value of 
log(N/H) + 12 = 8.30 + 0.30 dex from our analysis in Sect.g] 
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Fig. 7. Comparison of the spectra of GHV-62024 (black) and NGC346-010 (red). 

reach the derived value of 8.3, although they are consistent with 
the lower half of the error bar. We have to take into account that 
the highest possible N enrichment is limited by the initial abun- 
dances of C and O (which decrease with the metallicity), and 
that the C and O abundances could scale differently from solar; 
alternative ly, there might still be some p roblems in our N m line- 
synthesis dRivero-Gonzalez et al. 11201 lah . which might add to 
the quoted error budget. The He abundance is also higher than 
predicted for a model with an initial mass Mq = 30 M Q at the 
chosen initial rotational velocity, but would be consistent with 
that of a model of 40 M Q , for which Yu e = 0.20 is predicted 
at the temperature obatined for GHV-62024, or with that of a 
slightly faster rotating model. 

Our results qual it itatively agree with those from 
iRivero-Gonzalez et all (1201 lbj) for the LMC, who also 
obtained very high N and He abundances for a large fraction 
of their LMC sample. GHV-62024 would very nicely fit in 
Rivero-Gonzalez et al. Fig. 8 (lower panel), where they show a 
clear correlation between the abundances of these two elements. 
This indicates that in both cases the N and He enrichment has 
the same origin, presumably the exposition at the surface of the 
products of the CNO bi-cycle taking place in the stellar interior. 

Using the same tracks, we can estimate the stellar age to be 
less than 4.5 Myr, in qualitative agreement with our expectations 
for a young region. In Fig. [10] we plot the behaviour of the ro- 
tational velocity versus the effective temperature in the models. 
Interestingly, we see that the rotational velocity is predicted to 
remain at high values until much later in the evolution of the 
star (the effective temperature acts here as a temporal variable). 
This is fully consistent with our expectation that the star is a fast 
rotator. Therefore, we think that it is possible to conclude that 
GHV-62024 shows strong N and He abundances that are com- 
patible with a high initial rotational velocity still present in the 
star, although other possibilities (such as an interactive binary 
history) cannot be discarded with the present data. 

However, if the star still retains its high rotational velocity 
(note in Fig. [10]that for the model with 30 M the rotational ve- 
locity increases during the early evolution), this could mean that 
2D effects may become important. In that case, our ID models 
would be too simplified and orientation effects could play a role. 
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Fig. 6. Wind-momentum luminosity r elationship for GHV- 
62024 (upper green triangle: value from iHerrero et alH (1201 ll) : 
middle green triangle: this work; lower green triangle: this work, 
considering a wind terminal vel ocity 30% lowe r than q uoted in 
Tab. [2} compared to data from Moki em et al. I (l2007bl) for the 
MW (red circles) and SMC (blue squares). The SMC stars cited 
in text were marked with an open circle. From left to right they 
are NGC346-010, AzV296, AzV8 3 and NGC346-001 . Th e solid 
lines represent the WLR fitted by Mo kiem et al. I d2007bl) with- 
out considering clumping for MW and SMC objects, and the 
black dashed line corresponds to the SMC relation shifted to 
Z = O.13Z . 



This strong nitrogen enhancement (and the similarly high He 
abundance, Yu e = 0.18 + 0.04) cannot be explained with non- 
rotating models, which only predict N and He enhancement 
much later, in the red supergiant phase. Therefore, we lo oked 
at the set of rotating models calculated bv lBrott et al. I (1201 lb for 
the SMC metallicity. Fig. [9] presents the predictions of models 
with y rot =410 km s compared with the value derived from 
the observations. The models predict a slightly lower maximum 
N abundance (except for homogeneous evolution) and never 
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Fig. 8. Position of GH V-62024 in the HR D iagram. Evolutionary 
tracks are taken from iBrott et a"L~l (1201 lb . with an initial rota- 
tional velocity of 410 km s . 
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Fig. 9. Variation of the abundance of N along the stellar life- 
time (as given by the evolution in temperature). Represented in 
the ordinate axis is \og{ u NIH) + 12. Models with a zero ini- 
tial rotational velocity only change their nitrogen abundance in 
the red supergiant phase, whereas models with an initial rota- 
tional velocity of 410 km s _1 suffer an early nitrogen enhance- 
ment. Evolutionary models have masses of 60, 50, 40, 35, 30, 20, 
15 and 12 M , from left to right. The position of GHV-62024 is 
represented by the dashed box. 



A proper comparison would have be made between 2D evolu- 
tionary and atmospheric models, which is a work for the future. 
However, we are confident that even if the high rotational ve- 
locity is confirmed, the analysis and comparisons described here 
are a good representation of the average properties of the star. 

7. Conclusions 

We have analysed GHV-62024, an 06.5 Illf star in the low- 

metallicity galaxy IC 1613. 

Our preliminary analysis of this star dHerrero et al. I SoTTb 
resulted in a too high MWM rate for the star, typical of Milky 
Way metallicity stars. This would have posed a problem for the 
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Fig. 10. Evolution of the rotational velocity for evolutionary 
models with initial masses of 30 (black), 35 (blue, dot-dashed 
line) and 20 (red , dashed line) solar masses, according to 
Brot t et al. I (1201 ll) . The temperature range derived for GHV- 
62024 is indicated by the hatched region. 



theory of radiatively driven winds. However, in our new analysis 
presented in this work we conclude that the MWM is consistent 
with a lower-than-SMC metallicity, but we derive a f3 exponent 
of 2.0 for the wind velocity law. This is higher than the usual 
value for stars of th is spectral type. Following a suggestion by 
iHillier et al. I (|2003), we hypothesize that the star is a fast rotator 
seen from low inclination angle and that the high ft would then 
be an artefact caused by the superposition of light emitted from 
regions with different wind properties. 

When compa r ing ste llar evolutionary models with rotation 
from Bro tt et al. I (1201 ll) we see that the star shows the well- 
known mass-discrepancy. We derive high N (log(N/H)+12.= 
8.30+0.30) and He (Y= 0.18+0.04) abundances, which can 
qualitatively be explained by a strong mixing with processed 
CNO products. The N abundance i s higher than the maximum 
value predicted by the tracks from IBrott et al. I (1201 ll) for non- 
homogeneous evolution, although it is consistent within the er- 
rors, and the He abundance is consistent with the model predic- 
tion either for a higher mass than the initial stellar mass derived 
from the position of the star in the HRD or for a faster initial ro- 
tational velocity. However, with an initial rotation of 410 km s _I 
the models predict that the star can reach high abundances while 
still keeping the high rotational velocity. This is therefore consis- 
tent with the above hypothesis that GHV-62024 is a fast rotator 
seen from a low inclination angle, whose velocity field proper- 
ties would have been modified towards high ft values without 
modi fying the average W LR (in agreement with lMaederll2001l 
and lPetrenz & Puis Il2000t) . 

However, a high rotational velocity could give rise to 2D ef- 
fects, not included in the atmosphere or evolutionary models. If 
the high rotational velocity is confirmed, the analysis and com- 
parisons carried out in this work should be considered as a first 
step towards a more complete, multidimensional analysis. 
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